A CHARACTERIZATION OF ITERATIVE EQUATIONS BY 
THEIR COEFFICIENTS 
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Abstract. An expression for the coefficients of a linear iterative equa- 
tion in terms of the parameters of the source equation is given both 
for equations in standard form and for equations in reduced normal 
form. The operator generating an iterative equation of a general order 
in reduced normal form is also obtained and some other properties of 
iterative equations are established. In particular, a simple necessary 
and sufficient condition for an equation to be iterative is given for the 
general fourth-order linear equation solely in terms of its coefficients. 



1. Introduction 

It is well-known [lj that linear ordinary differential equations (lodes) 
of order one or two can all be reduced by a local diffeomorphism of the 
(x, y)-plane to the canonical form y' = and y" = 0, respectively, and 
that this is not the case for equations of a general order n > 2. Lie [2] 
showed that a differential equation of a general order n > 2 is equivalent 
(by a local diffeomorphism of the plane) to the equation y^ = 0, which we 
shall henceforth refer to as the canonical form of the linear equation, only if 
its symmetry algebra has the maximal dimension n + 4. In a much recent 
paper, Krausse and Michel [3] proved the converse of this statement and 
also showed that a LODE of order n > 2 has a symmetry algebra of maximal 
dimension if and only if it is iterative. Linear iterative equations are the 
iterations ty n y = of a linear first order equation, of the form 

^y = r{x)y' + s(x)y = 0, $ n y = V^^Wy. (1.1) 

In these iterations, the equation r(x)y' + s(x)y = is termed the source 
equation, while the functions r(x) and s(x) are referred to as its parameters. 
Despite the unique symmetry properties of these iterative equations, many 
of their properties are still not known, and Mahomed [lj gave a listing of 
these equations for the orders three to five. We extend this list to equations 
of a general order, in both the standard form and the associated reduced 
normal form, and determine the operator generating the linear iterative 
equation of any given order in reduced normal form. Other properties of 
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iterative equations are also established and in particular a simple criterion 
for a fourth-order linear equation to be iterative is given solely in terms of 
its coefficients. 

2. Equations in the general linear form 

By replacing the dependent variable y = y(x) by y + y p , where y p is a 
particular solution of the inhomogeneous equation, we may assume without 
loss of generality that a linear iterative equation of a general order n has the 
form 

= K° n yW + K\ y^ + K 2 n y^ + ■■■+ K n n ~ l y> + K n n y = K. (2.1) 

2.1. Case where \& = ^ + s. For the sake of clarity we first consider the 
case where the differential operator ^ = rd/dx + s is much simpler, with 
r ee r{x) = 1. It is clear that in this case, the operator \& leaves invariant the 
leading coefficient, and thus we have K® = 1 in this case. Using the formula 
for ^f n y in (jl.ip gives the recurrence relations 

K l n = K l n _ l + s. (2.2a) 

K = K-i + i {kS) + = K-i + *K~l (2.2b) 



dx 

(for j = 2, . . . , n — 1 and n > 2). 

K = Y X ( K n~D + Si ^-1 = * K nZl (2.2C) 

Setting 

-fC^ = 0, for j < or j > to, and iT^ = 1, for m = j = (2.3) 
reduces the recurrence equations (|2.2[) to the single equation 

Ki = Kl_ 1 + ^Kiz\, 0<j<n, Vn>l. (2.4) 
Solving the recurrence relations (|2.4p together with the initial conditions 

= y' + sy = y' + K\y 
^ 2 y ee y" + 2sy' + (s 2 + s')y = y" + K\y' + K 2 y 

readily gives 

K = * n_1 s> 

for all n > 1. Using these two equalities, and setting 

f = 1, for every function / = f(x), (2-5) 
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one readily sees by induction on j and n that 

* = to=*...„.»d„>i. 

We have thus obtained the following result for the case where r = 1. 

Theorem 1. When the operator ^> generating the iterative equation has the 
form = d/dx + s, that is, when it depends on the single function s, the 
coefficients Kh of the iterative equation of a general order n are given by 



Kl 



Vj = 1, 



(2.6) 



and the iterative equation of a general order n is therefore given by 



v n y = y 



(n) 



3=1 



n 



(2.7) 



Note however that (|2.6p is also valid for j = 0. 

2.2. Case where the source equation depends on both parameters. 

In this case we have = rd/dx + s, where the parameters r, and s are given 
functions, and this is the most general case. Using the definition for $> n y 
given in (|2.ip as well as the conventions for Kh set in (12. 3j) show that the 
coefficients Kn of the iterative equation (|2.1[) satisfy the recurrence relations 



K 3 n = + ^K J n _\, for < j < n, and n > 1, 



(2i 



which naturally reduce to (12. 4|) for r = 1. Setting j = or j = n in (J2T 
readily gives by induction on n the identities 



n-l. 



for all n > 1. 



(2.9) 



Applying (|2,8p recursively and using the conventions set in (|2.3p give a new 
recurrence relation 



K i = T, rn ~ k * K £ 1 v fori = 0,. 
fc=j 



, n and n > 1 . 



(2.10) 



Although Eq. (I2.10p does not provide the required expression for Kh in 
terms of the parameters r and s, it represents an algorithm for the compu- 
tation of the coefficients Kh for all possible values of n and j. For instance, 
using (I2.10|) with j = 1,2 readily gives 



Ki 



Kt 



„n-2 



ns + 



n \ T / n 
2* S+ 3 



3sr' + rr" + 



3n — 5 



J'2 



(2.11a) 
(2.11b) 
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To obtain the general expression for K n using a recurrence relation relating 
them, we rewrite (|2.10p in the form 

n 

Ki=T, r(k - j) * K i~l +j -r ( 2 - 12 ) 

k=j 

Then, using ()2.3p and (12. 9h . the following formulas are successively obtained: 



Ki 



J2r 




k=l 




n 


n— k 2 + l 


E 


E 


k 2 =2 


fci=i 


n 

E 


n— fc 3 +5 

E 


fc 3 =3 


k 2 =2 




r r 



fcl=l 



(2.13a) 
(2.13b) 

(2.13c) 



Continuing this process with two more iterations by computing and K\ , 
a clear pattern for the general coefficient Kh emerges, and to write down 
this expression we introduce some notations. For n > 1 and < % < j < n, 
set 



u=i+l 

Mi = MiO') =n + 
0' 



where k u £ 1* 
j 



ay = n + x } 



Poj = M . 



(2.14a) 

(2.14b) 
(2.14c) 



We have the following result for the general case, where ^ = rd/dx + s. 

Theorem 2. In terms of the parameters r and s of the source equation, the 
general coefficient K n of the iterative equation (12. ip has the form 



M 2 Mi 

*i = E E "E E 

fcj=j kj—i=j—l k 2 =2 fci = l 



(2.15) 



/or n > 1, and 1 < j < n, and where the expressions for {3ij,Mi, and ay are 
gu>en 6?/ f|2.14j) . 



For the sake of clarity it would be useful to verify explicitly that (|2.15p 
reduces indeed to (12.61) for r = 1. 



Proposition 1. Eq. H2.15H reduces as expected to (I2.6f) for r = 1. 

Proof. When r = 1, the general term in the summation (I2.15D clearly reduces 
to \|W ■ 1 = & 3 s, and since this expression does not depend on the running 
indices ki, kz, . . . kj to prove the proposition, it suffices to show that the 
total number 

Mj Mj_i M 2 Mi 

^ = E E •••EE 1 ( 2 - 16 ) 

kj=j fcj_i=J — 1 ^2=2 fci=l 

of terms in this summation is precisely (™). For j = 1, and j = 2, it clearly 

follows from (j2. 13aH and (I2.13bj) respectively that P n ^\ = (™) , and P n ^i = (2) ■ 
It also follows from (|2,9p that P n $ and P n ^ n also satisfy the required property 
for n > 1. Let n>2 and assume that P v ,j-i = (j-i)> f° r 1 < w < n and 
0<j — l<f<n. Then it follows from (|2. 12|) and the linearity of ^ that 
P n j = Ylk=j { n ~ l j-{~ 1 ) ■ Setting n — k + j — 1 = q and j — 1 = m gives 

= E (™) = ( ( m + 1 ) = l) ' 

q=m x 7 x / \" / 

and this completes the proof by induction of the required property for P n ,j- 

'□ 

Although Proposition Q] gives a verification of the validity of the compli- 
cated formulas (|2.15p at least for the simpler case r = 1, these formulas can 
be slightly simplified, by a suitable change of variables. Indeed, if in (|2,15p 

we set 

(k j -j = n-P j , for j>l 

\ki-i = P i+ i - Pi - 1, for i = 1, . . . , j - 1 

This reduces the expression for Mj in (|2.14p to 

Mi = P i+1 - 1 (2.18) 

In particular ctj = M Q is reduced to 

Oj = Pi - 1 (2.19) 

Since Mj = n for all j, thanks to (I2.18p . we may rewrite (I2.17P as 

k i -i = M l -P i , for i = l,...,j. (2.20) 

Consequently, thanks to (|2.18p . (I2.19P and (I2.20p . and after renaming the 
Pi, we may rewrite (I2.15P in the slightly simplified form 

n Ay — 1 £,'2 — 1 

k> y, E ••■EE 

k 2 =2 fci=i (2-21) 
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3. Equations in reduced normal form 

By diving through the general n-th order linear iterative equation ty n y in 
(|2.1j) by = r n , it can be put in the form 

y (n) + B l y(n-i) + ... + BI y^-n + ■ ■ ■ + Bl y = 0, (3.1) 

where B n = K n /r n . This is the standard form of the general linear iterative 
equation with leading coefficient one. It is well-known that (|3.ip can be 
transformed to the normal form (in which coefficient of y( n_1 ) has vanished) 
by a change of the dependent variable of the form 

y^yexpQ^ B l n (v)dv^j . (3.2) 

However, this amounts to the requirement that B n = 0, i.e. that K n = 0. 
Therefore, an n-th order linear equation in reduced normal form is iterative 
if and only if it has the form 



,(n) + A l y{n-2) + ... +A l y (n-j) + ... A n n y = (3.3a) 



where 



K 3 

A3 - ttH. 

A n - n 



(2 < j < n), (3.3b) 

Kl=o 



and where K 3 n is given by (f2TT5]) . It follows from (|2.11ap that setting K\ = 
amounts to setting 

s = -~(n-l)r', (3.4) 

and this shows why any iterative equation in normal form can be expressed 
in terms of the parameter r alone. Moreover, thanks to (|3.3b|) the coefficients 
A J n inherit all of the characterization (I2.15P obtained for the coefficients K n 
of iterative equations in standard form, and the corresponding expression 
for the An has up to the factor l/r n , exactly the same form as that for the 
K n because the parameter s does not appear explicitly in Eq. (|2.15p . In 
addition, using the same algorithm (e.g. (|2.1U|) ) obtained for the K n , one 
can readily compute An for all n > 2 and 2 < j < n. For instance, using the 
expression for K\ in (|2.11bp together with (|3.3bj) . one readily sees that 

^_fn + l\ Af , a2 r' 2 -2rr" 



K=[ 3 )A(r), where A{r)=A£ = — 2 . (3.5) 



Theorem 3. Let the differential operator & n be given by 



(3.6) 



Kl=0 
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Then the equation & n y = is exactly (I3.3ap . that is <3? n generates the (most 
general) linear iterative equation of an arbitrary order n in normal form. 

Proof. Indeed, *$> n y generates the linear iterative equation of general order 
n in standard form and (l/r n )^ n generates the same equation with leading 
coefficient 1, while setting K\ = corresponds as already noted to reduce 
the latter equation to its normal form. □ 



4. Applications 

Although formula f|2. 15j) in Theorem [2] gives a characterization of the 
coefficients of a linear iterative equation of a general order in terms of the 
parameters of the source equation, and therefore a characterization of the 
linear iterative equation itself, in practice a linear equation is given solely 
in terms of its coefficients, and without reference to any source equation or 
the parameters thereof. We thus need a characterization of these iterative 
equations that relies solely on the coefficients of the equation. Using the 
operator <5 n of Theorem [3] we can easily generate an iterative equation of 
a general order n, and we let the generated equation be in the form (|3.3ap . 
By analyzing the coefficients of the latter equation we readily see, at least 
for low order equations, that they all depend only on the coefficient A\ and 
its derivatives. For instance, for n = 3 or 4, if we set A\ = 03 and A\ = a^, 
then iterative equations of orders 3 and 4 take on respectively the forms 



V'" + a 3 y' + -a' 3 y = (4.1) 
yW + a A y" + a' A y> + (±aT A + ^af) y = (4.2) 

Note that conversely, any equation of the form (|4.ip or (|4.2p is iterative, be- 
cause we can always solve for r the equation A\ = ( n J }A(r) appearing in 
(|3.5p . together with (|3.4p to find the parameters of the corresponding source 
equations. Thus each of the equations (|4.ip and (|4.2p characterizes the nor- 
mal form of the iterative equation of the corresponding order. However, the 
most general characterization is to be given for equations in standard form 
(|3.ip . and the most practical characterization should be expressed explicitly 
in terms of the coefficients of the equation. 

Let a third-order linear lode be given in the form 

y'" + c 2 y" + c iy ' + c y = 0, (4.3) 

where the coefficients Cj for j = 0,1,2 are all functions of x. We may 
thus assume that its reduced normal form is given by (|4.ip . Reverting back 
the latter reduced equation to the corresponding equation in standard form 
using the inverse of the transformation of type (13. 2j) , and dropping the 
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subscripts in the expression of the coefficients in (14. ip gives 

w'" + c 2 w" + (a + c' 2 + c|/3) w' 

1 ' (4-4) 

+ — (27a' + 18c' 2 ' + 18ac 2 + 18c' 2 c 2 + 2c|) iu = 0, 

where w is the new depend variable. Letting w = y in (|4.4p and then 
equating its coefficients with those of (14.3P shows that 



a = ci — (c' 2 + c 2 /3) (4.5a) 
Co = [27a' + 18c' 2 ' + 18ac 2 + 18c' 2 c 2 + 2cf] . (4.5b) 

Using (I4.5ap for the expression of a and its derivatives and substituting the 
result in (I4.5bp gives 



54c - 18cic 2 + 4eij ~ 2 7ci + 18c 2 c' 2 + 9c 2 = 0. (4.6) 

The latter equation is the well-known characterization due to Laguerre [5] 
and Lie [B] of linear third order equations that can be reduced to the canoni- 
cal form y'" = by a point transformation. However, the methods they used 
to derive this conditions were different from the one used here, and which is 
based on the characterization of iterative equations given by fl2. 151) and the 
generating operator <3? n in (]3.6p . Clearly, due to the result already cited of 
Krause and Michel [3] relating iterative equations and equations reducible 
to canonical form by point transformations, (14. 6p is also a characterization 
of linear third order equations which are iterative. 

We now move on to consider the case of a fourth order linear equation 
given in the form 

2/ {4) + c 3 y"' + c 2 y" + c iy ' + c y = 0, (4.7) 

and we are interested in deriving necessary and sufficient conditions on the 
coefficients Cj (j = 0, . . . , 3) for the equation to be iterative. Here again, we 
may assume that the reduced normal form of the equation is given by (14.21) . 
Transforming back the latter equation to its standard form, dropping the 
subscript in and equating the coefficients of the resulting equation with 
those of fl4.TR gives 

« = ^(-3ci + 8c 2 -12c' 3 ) 

ci = ^ac 3 + —4 + a> + + 4 

6400c = 576a 2 + 400a (4cg + c§) + 80 (l5cf + 24a" + 204") 
+ 1600 (a + 4) C3 + 600</ 3 c§ + 25^. 



(4.8a) 
(4.8b) 

(4.8c) 
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Substituting in (|4.8bj) and (|4.8cj) the expressions for a and its derivatives 
given by fl4.8a[) yields 

= 4c 2 c 3 -4 + 8c 2 - 6c 3 c 3 - 4c 3 ' - 8ci (4.9a) 

= 1600c - 144c 2 , + llcj - 400c 3 c' 2 + 288cfcg + 336C3 2 

(4 9b) 

+ 8c 2 (4 + 4 4) - 480c 2 + 56OC3C3 + 320c 3 ". 
We thus have the following result. 

Theorem 4. A Linear fourth order equation of the general form (j4.7|) is 
iterative if and only if its coefficients satisfy the system of two equations 

na. 



Proof. As the system (|4.9p was obtained under the assumption that (|4,7p 
is iterative, we only need to prove conversely that the equation is iterative 
whenever its coefficients satisfy (|4.9p . The reduced normal form of (14. 7p 
has, after the substitution of the expressions for cq and c\ given by (|4.9j) in 
terms of 02,03, and their derivatives, the form 

w {A) + Q 2 w" + Q x w' + Q w = (4.10a) 
where 

Q 2 = c 2 -^(c| + 4c' 3 ) (4.10b) 

Qi = 4 - ^(c 3 c 3 + 24) (4.10c) 



+ ^( 27c t + 216c !4 + 640c 2 " 480c 3 c 3 ' - 960c 3 "). 
The coefficients Qj thus obtained clearly satisfy the conditions 
Qi = Q' 2 and Q = (Aq' 2 ' + Aq2 } 



(4.10d) 



prescribed by (|4.2p . and this completes the proof of the theorem. □ 

5. Concluding remarks 

Some of the results that we've obtained in this paper include the expres- 
sion given in (|2.15p for the coefficients of an iterative equation of a general 
order in standard form (|2.ip in terms of the parameters of the source equa- 
tion, and the similar expression in (j3.3b|) for iterative equations in reduced 
normal form (|3.3ap . We have provided a simple algorithm (I2.10P for the 
calculation of these coefficients for given parameters of the source equation, 
and the operator $ n obtained in Theorem [3] generates the iterative equation 
of an arbitrary order n in normal form. Given the already cited result of 
Krause and Michel [3j according to which an equation is iterative if and only 
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if it is reducible by a point transformation to the canonical form, the char- 
acterization of iterative equation we have given in (14.9H is also an extension 
to the order four of the characterization (I4.6P due to Laguerre [5] and Lie 
|6] of linear third-order equations that are reducible to the canonical form. 
But unlike the characterization (I4.6P which consists of a single equation, the 
characterization (14. 9f) for fourth-order equations consists of a system of two 
equations, and from the process of their derivation it should be expected 
that for equations of order n the corresponding characterization will consist 
of a system of n — 2 equations. 
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